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Abstract. In the context of supersymmetric models in which small Dirac neutrino masses 
are generated by supersymmetry breaking, a mainly right-handed (RH) mixed sneutrino 
can be an excellent cold dark matter (DM) candidate. We perform a global analysis of 
the Minimal Supersymmetric Standard Model (MSSM)+RH neutrino parameter space by 
means of Markov Chain Monte Carlo sampling. We include all relevant constraints from 
collider and dark matter searches, paying particular attention to nuclear and astrophysical 
uncertainties. Two distinct cases can satisfy all constraints: heavy sneutrino DM with mass 
of order 100 GeV, as well as light sneutrino DM with mass of about 3-6 GeV. We discuss 
the implications for direct and indirect dark matter searches, as well as for SUSY and Higgs 
searches at the LHC for both, the light and the heavy sneutrino dark matter case. The light 
sneutrino case is excluded by the 125-126 GeV Higgs signal. 
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1 Introduction 

The nature of dark matter (DM) [1—3], the nature of the physics stabilizing the electroweak 
scale [4-6], and the origin of neutrino masses [7-9] range among the most challenging open 
problems in particle physics. Huge efforts on both the experimental and theoretical sides 
are undertaken worldwide to shed light on these questions and eventually unravel a more 
fundamental theory beyond the current Standard Model (SM) of electroweak and strong 
interactions. 

Weak-scale supersymmetry (SUSY) [10-13] is a fascinating option for such a theory 
beyond the SM, potentially addressing all three of the above open problems. Indeed SUSY, 
if realized at the TeV scale, elegantly solves the gauge hierarchy problem and, if R-parity is 
conserved, provides an excellent particle dark matter candidate: the lightest supersymmetric 
particle (LSP). Moreover, in a certain class of models, small neutrino masses may naturally 
arise from F-term SUSY breaking [14, 15]. In addition to providing an explanation for 
neutrino masses, this class of SUSY models offers a particular DM candidate: a mainly 
right-handed (RH) mixed sneutrino. 
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Mixed sneutrinos as thermal DM are indeed a very interesting alternative to the con- 
ventional neutralino LSP of the Minimal Supersymmetric Standard Model (MSSM). They 
have received much attention recently, in part because of their intriguing phenomenology and 
in part because they provide a possibility for light SUSY DM below 10 GeV. Many studies 
of sneutrino DM have been performed in models with extra singlets, extra gauge groups or 
models with Majorana neutrino masses, see e.g., [16-39]. In our work, we concentrate instead 
on the MSSM+RH neutrino model [14, 15] with only Dirac masses for neutrinos. The phe- 
nomenology of this model was investigated in detail in [14, 40]. Indirect detection signatures 
were discussed in [41, 42], implications for O^/Odm in [43], and LHC signatures in [44, 45]. 

The crucial point of this model is that one can have a weak-scale trilinear sneutrino 
coupling Ay that is not suppressed by a small Dirac-neutrino Yukawa coupling. It can 
hence induce a large mixing between left-handed and right-handed sneutrinos even though 
the Yukawa couplings may be extremely small. The lightest sneutrino can thus become the 
LSP and a viable thermal DM candidate. Note that the mainly RH sneutrino LSP is not 
sterile but couples to SM gauge and Higgs bosons through the mixing with its LH partner. 
Sufficient mixing provides efficient annihilation so that the sneutrino relic density Qh 2 ~ 0.11 
as extracted from cosmological observations [46]. 

Direct detection (DD) experiments however pose severe constraints on Dirac or complex 
scalar, i.e. not self-conjugated, DM particles because the spin-independent elastic scattering 
cross section receives an important contribution from Z exchange, which typically exceeds 
experimental bounds. In the mixed sneutrino model, this cross section is suppressed by the 
sneutrino mixing angle. Therefore, on the one hand a viable sneutrino DM candidate requires 
enough mixing to provide sufficient pair-annihilation, on the other hand the mixing should 
not be too large in order not to exceed the DD limits. 

In [40] , some of us explored the parameter space of the Dirac sneutrino DM model where 
these conditions are satisfied for light sneutrinos with a mass below 10 GeV. This mass range 
was motivated by hints of DM signals in DD experiments [47, 48]. In the present study, we 
explore a much wider range of masses, considering both light DM below 10 GeV as well as 
heavier DM of the order of 100 GeV. Moreover, we explore the parameter space by means 
of Markov Chain Monte Carlo (MCMC) sampling, using Bayesian statistics to confront the 
model predictions with the data. In taking into account the limits from DD experiments 
we pay special attention to uncertainties stemming from astrophysical parameters (local 
DM density and velocity distribution) and to uncertainties in the quark contents of the 
nucleons (relevant in particular when there is a large Higgs-exchange contribution). Finally, 
we consider the impact of the LHC results, which push the masses of squarks and gluinos 
above the TeV scale. Our results are presented as posterior probability densities of parameters 
and derived quantities, in particular of the DM mass and direct and indirect detection cross 
sections. 

The paper is organized as follows. In Section 2, we briefly recall the main features of 
the mixed sneutrino model. In Section 3, we then describe in detail the setup of our MCMC 
analysis. The results of this analysis for light sneutrino DM are presented in Section 4.1, and 
for heavy sneutrino DM in Section 4.2. Conclusions are given in Section 5. Appendix A 
contains a discussion of the uncertainty in the effective degrees of freedom in the early 
Universe, Appendix B contains summary tables for our results, and Appendix C illustrates 
the prior dependence by comparing to log-prior results. 
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2 Framework 



The framework for our study is the model of [14, 15] with only Dirac masses for neutrinos. 
In this case, the usual MSSM soft-breaking terms are extended by 

A£ soft = m%]Ni\ 2 + A^LiNiHu + h.c. , (2.1) 

where m 2 ~ and A u are weak-scale soft terms, which we assume to be flavor-diagonal. Note 
that the lepton-number violating bilinear term, which appears in case of Majorana neutrino 
masses, is absent. Neglecting the tiny Dirac masses, the 2x2 sneutrino mass matrix for one 
generation is given by 

f m% + \rn? 7 cos2f3 ^sA v v sin /A 
ml = L i \ . a V2 2 • ( 2 -2) 

Here m~ is the SU(2) slepton soft term, v 2 = v 2 + v\ = (246 GeV) 2 with v\ 2 the Higgs 

I j 

vacuum expectation values, and tan/3 = Vilv\. The main feature of this model is that m~, 

Li 

m 2 ~ and Ay are all of the order of the weak scale, and Ap does not suffer any suppression 
from Yukawa couplings. In the following, we will always assume < mj^ so that the lighter 
mass eigenstate, ui, is mostly a ur. This is in fact well motivated from renormalization group 
evolution, since for the gauge-singlet N the running at 1-loop is driven exclusively by the Ay 
term: 

dm 2 - a 

^f = T^ < 2 ' 3 > 

while 

2 



dm- 



Li _ ft.ccM i._„\ , 2 2 



. , - (MSSM terms) + — ^ A£ . (2.4) 
at 1d7t z 1 



The renormalization group equation (RGE) for the ^4-term is: 

Here, g\ and g2 are the U(l) and SU(2) gauge couplings, and y% and are the top and 
charged lepton Yukawa couplings. 

A large A v term in the sneutrino mass matrix will induce a significant mixing between the 
RH and LH states, 

\V2 ) \ sm V v cos 8y J \ vl J m z ~ 2 - mr~ 

and a sizable splitting between the two mass eigenstates v\ and V2 (with m Ul < m U2 ). 

One immediate consequence of this mixing is that the mainly RH state, i>\, is no longer 
sterile. However, its left-handed couplings are suppressed by sin By . This allows the v\ to 
have a large enough pair-annihilation rate to be a viable candidate for thermal dark matter, 
while at the same time evading the limits from direct dark matter searches [14, 40, 44]. A 
mainly RH v\ as the LSP will also have a significant impact on collider phenomenology, as it 
alters the particle decay chains as compared to the "conventional" MSSM. Moreover, it can 
have a significant impact on Higgs phenomenology: first, a light mixed sneutrino can give a 
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large negative loop correction to m h o which is oc \A U \ 4 [40]; second, a large A v can lead to 
dominantly invisible Higgs decays if m Ul < m h o/2. 

In the following, we will assume that electron and muon sneutrinos are mass-degenerate, 
m Uie = mp i with %= 1,2. Moreover, by default we will assume that the tau-sneutrino, v\ T 
is lighter than the D\ e and is the LSP. This is motivated by the contribution in the running 
of the j4-term coming from the Yukawa coupling, see eq. (2.5). In this case, we take m^, 
mc 2 , sin#;> and tan/3 as input parameters in the sneutrino sector, from which we compute 
mg, rrifr, Ay (all parameters are taken at the electroweak scale). 

3 Analysis 
3.1 Method 

We choose to confront the sneutrino DM model to experimental constraints by means of 
Bayesian inference. In this kind of analysis, one starts with an a priori probability density 
function (prior PDF) p(6\A4) for the parameters 6 = {#i... n } of the model Ai, and some 
experimental information enclosed in a likelihood function p(d\6,M) = C{9). The purpose 
is to combine these two pieces of knowledge, to obtain the so-called posterior PDF, possi- 
bly marginalized to some subset of parameters. Splitting the parameter set as 9 = (ip,A), 
Bayesian statistics tells us that the posterior PDF of the parameter subset ip is 

p(il>\M) oc J d\p(iP,X\M)C(ip,X). (3.1) 

That means one simply integrates over unwanted parameters to obtain the marginalized 
posterior PDFs. These unwanted parameters can be model parameters, but can also be 
nuisance parameters. 

In this work, we evaluate posterior PDFs by means of a Markov Chain Monte Carlo 
(MCMC) method. The basic idea of a MCMC is setting a random walk in the parameter space 
such that the density of points tends to reproduce the posterior PDF. Any marginalisation 
is then reduced to a summation over the points of the Markov chain. We refer to [49, 50] 
for details on MCMCs and Bayesian inference. Our MCMC method uses the Metropolis- 
Hastings algorithm with a symmetric, Gaussian proposal function, basically following the 
procedure explained in [51]. We use uniform (linear) priors for all parameters. The impact of 
logarithmic priors in the sneutrino sector is presented in Appendix C. For each of the scenarios 
which we study, we run 8 chains with 10 6 iterations each, and we check their convergence 
using the Gelman and Rubin test with multiple chains [52], requiring \fl\ < 1.05 for each 
parameter. First iterations are discarded (burn-in), until a point with log(£) > —5 is found. 

The likelihood function C can be constructed as the product of the likelihoods asso- 
ciated to the N observables 0«, 

N 

C = Y[Ci. (3.2) 

i=l 

Available experimental data fall into two categories: measurements of a central value, and 
upper/lower limits. In the former case, the central value O exp comes with an uncertainty 
given at some confidence level CL. It is reasonable to assume that the likelihood function for 
this kind of measurement is a Gaussian distribution, 

Ct = M{0 - O cxp , AO) = exp ( ~ { ° m °y )2 ) ■ (3-3) 
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Here AO is the uncertainty at la. For combining experimental and theoretical uncertainties, 
we add them in quadrature. When O exp is a (one-sided) limit at a given CL, it is less 
straightforward to account for the experimental uncertainty. Taking a pragmatic approach, 
we approximate the likelihood by a smoothed step function centered at the 95% CL limit 

O ex p, 95%) 

A = F(0,O exPi 95%) = — 7-—— — ^ waTSi ' ( 3 - 4 ) 

1 + exp[±(0 - O exPj95% )/AOJ 

with AO = 1% x O exp 95%. The ± sign in the exponent is chosen depending on whether we 
are dealing with an upper or lower bound: for an upper bound the plus sign applies, for a 
lower bound the minus sign. Using a smeared step function rather than a hard cut also helps 
the MCMC to converge. 

Finally, when the \ 2 °f the limit is available (this will be the case for the direct detection 
limits), we compute the likelihood as Li = e~ Xi ^ 2 . 

To carry out the computations, we make use of a number of public tools. In particular, 
we use micrOMEGAs 2.6. c [53, 54] for the calculation of the relic density and for direct and 
indirect detection cross sections. This is linked to an appropriately modified [40] version of 
SuSpect 2.4 [55] for the calculation of the sparticle (and Higgs) spectrum. Decays of the 
Higgs bosons are computed using a modified version of HDECAY 4.40 [56], and Higgs mass 
limits are evaluated with HiggsBounds 3.6.1beta [57, 58]. Regarding the computation of 
the direct detection limits, we make use of a private code described in section 3.4.2. 

3.2 Parameters of the model 

We parametrize the model with 12 parameters as follows. The sneutrino sector is fixed 
by three parameters per generation (the two mass eigenvalues m^, and the mixing 
angle sin 9^, or the soft breaking parameters m^, m^, Ay) plus tan/3. Assuming degeneracy 
between electron and muon sneutrinos, this gives seven parameters to scan over. The soft 
term for the LH sneutrino, m^, also defines the mass of the LH charged slepton (of each 
generation); the remaining free parameter in the slepton sector is m^, the soft mass of the 
RH charged slepton, which we fix by = toj for simplicity. 

The chargino-neutralino sector is described by the gaugino mass parameters M\ , M2 and 
the higgsino mass parameter [i. Moreover, we need the gluino soft mass M3. Motivated by 
gauge coupling unification, we assume [approximate] GUT relations for the gaugino masses, 
M3 = 3M2 = 6M1, so we have M2 and fi as two additional parameters in the scan. For 
stops/sbottoms we assume a common mass parameter 771,03 = m Q 3 = m u 3 = m 5 3 ' wm ch 
we allow to vary together with At (other trilinear couplings are neglected). The masses of 
the 1st and 2nd generation squarks, on the other hand, are fixed at 2 TeV without loss of 
generality. Finally, we need the pseudoscalar Higgs mass Ma to fix the Higgs sector. The 
model parameters and their allowed ranges are summarized in Table 1. 

The requirement of having enough sneutrino annihilation to achieve $7/i 2 ~ 0.11 while 
having a low enough scattering cross section off protons and neutrons to pass the DD limits, 
together with the constraints from the Z invisible width, splits the parameter space into two 
disconnected regions with sneutrinos lighter or heavier than Mz/2 (or more precisely, as we 
will see, < 7 GeV and > 50 GeV). We call this the "light" and "heavy" cases in the 
following. 

1 This assumption is central when applying the gluino mass limits from LHC searches. 
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i 


Parameter 


Scan bounds 




Pi 


light sneutrinos 


HND sneutrinos 


1 


m vlT 


[1, M z /2] 


[M z /2, 1000] 


2 




[m DlT + 1, 3000] 


[m 5lT + 1, 3000] 


3 


sin 6y T 


[0, i] 


[0, i] 


4 


fie 


[m 5lT +l, M z /2] 


K+l, 3000] 


5 




[m^+ 1, 3000] 


[m^+ 1, 3000] 


6 


sin 9p e = sin 9^ 


[0, i] 


[0, i] 


7 


tan j3 


[3, 65] 


8 




[-3000, 3000] 


9 


M 2 = 2Mi = M 3 /3 


[30, 1000] 


10 


m Q 3 = m % = m 5 3 


[100, 3000] 


11 


A t 


[-8000, 8000] 


12 


M A 


[30, 3000] 



Table 1. Parameters and scan ranges for the light and the heavy non-democratic (HND) sneutrino 
cases. All masses and the A-term are given in GeV units. In the heavy democratic (HD) case, the 
same bounds as in the HND case are applied for quantities i = 1-3 and 7-12, but entries 4-6 are 
computed from G ± 5%, € ± 5%, and Ap e g Ap T ± 5%, with a flat distribution, 
see text. 



In the "light" case, we assume that the r-sneutrino is the LSP, but the e/fi sneutrinos 
are not too different in mass from the r-sneutrino. More specifically, we assume that mo le lies 
within [m„ lT + 1 GeV, Mz/2], i.e. the tau sneutrino is the LSP and all the three sneutrinos 
are potentially in the region sensitive to the constraint on the invisible decays of the Z boson. 
The 1 GeV minimal mass splitting is a quite natural assumption considering the sensitivity of 
m^x t° small variations in A u , and suppresses co-annihilation effects (note that the degenerate 
case was previously studied in [40]). 2 

In the "heavy" case, we distinguish two different scenarios. First, in analogy to the 
light case, we assume that the r-sneutrino is the LSP, with m UlT G [Mz/2, 1000 GeV], and 
we allow m Ule to vary within [m UlT + 1, 3000] GeV. We call this the "heavy non-democratic" 
(HND) case in the following. Second, we also consider a "heavy democratic" (HD) case, in 
which m Ul , m„ 2 and sin 9 U of the 3rd and the lst/2nd generation are taken to be close to each 
other. As before, we use m UlT , m^ and sm9 Ur as input parameters, from which we compute 
m N ' m L an< ^ Av T - For the lst/2nd generation, we then take m^ G [m^ — 5%, m^ + 5%], 

G [mj — 5%, + 5%], and A Ve G [Ac T — 5%, A$ T + 5%] with a flat distribution. This 
way either D\ T or v\ e ^ can be the LSP; moreover v\ T and v\e,n can be almost degenerate. 
In the latter case, co-annihilations have a sizable effect. 3 Nevertheless it turns out that the 
results for the HND and HD setups are almost the same, so we will take the HND scenario 

2 We also performed MCMC sampling allowing mp lc > Mz/2 up to 3 TeV, keeping only the v\ T light, but 
the conclusions remain unchanged. So we will present our results only for the case mo lT < rnv 1( , < Mz/2. 

3 Note that if the electron/muon/tau sneutrinos are co-LSPs, this has important consequences for the relic 
density [40]. The e,fi,T sneutrino mass hierarchy moreover has important consequences for the LHC phe- 
nomenology (more electrons and muons instead of tau leptons from cascade decays), and for the annihilation 
channels for indirect detection signals. Furthermore, for a very light r-sneutrino, m^ lT < m T ~ 1.78 GeV, 
annihilation into a pair of tau leptons is kinematically forbidden, while for vie,^ of the same mass annihilations 
into electrons or muons would be allowed. 
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as our standard setup for the heavy case, see Table 1, and discuss only what is different in 
the HD case. 

3.3 Nuisance parameters 

Nuisance parameters are experimentally determined quantities which are not of immediate 
interest to the analysis but which induce a non-negligible uncertainty in the (model) param- 
eters which we want to infer. The Bayesian approach allows us to deal easily with nuisance 
parameters. In order to account for experimental uncertainties impacting the results, we 
choose 10 nuisance parameters, listed in Table 2. They fall into three categories: astro- 
physical parameters (related to dark matter searches), nuclear uncertainties (related to the 
computation of the DM-nucleon scattering cross section) and Standard Model uncertainties. 

In order to compute limits from direct detection experiments, we need to know the 
properties of the dark matter halo of our galaxy. We assume a Standard Halo Model, taking 
into account variations of the velocity distribution (vq, v esc ) and of the local dark matter 
density (/?dm)- To this end, we follow [59] and take the naive weighted average of the quoted 
values for each parameter (an alternative determination of pdm can be found in Ref. [60- 
62]). Note that considering vq and v esc as nuisance parameters is particularly important in 
the light DM case, because of its sensitivity to the tail of the velocity distribution; indeed a 
departure from the canonical value vq = 220 km/s may have a sizable impact on the direct 
detection limits at low masses. 

Turning to nuclear uncertainties, the Higgs exchange contribution to the elastic scat- 
tering cross section depends on the quark contents of the nucleons. The light quark contents 
can be determined via the ratio of the masses of the light quarks, m^/m^ and m s /m,i, and 
the light-quark sigma term a^N = {fn u + md){N\uu + dd\N)/2. Moreover, we need the 
strange quark content of the nucleon, a s = m s (N\ss\N) , which is actually the main source 
of uncertainty here. We take the latest results for a^N and a s from lattice QCD [63]. We 
stress that the new direct determinations of a s lead to a much lower value as compared to 
previous estimates based on octet baryon masses and SU(3) symmetry breaking effect. 

The Standard Model uncertainties that we include as nuisance parameters in the MCMC 
sampling are mt, the top pole mass, m;,(mj), the bottom mass at the scale nib i n the MS 
scheme, and a s (Mz), the strong coupling constant at the scale Mz- They impact the deriva- 
tion of the SUSY and Higgs spectrum. Moreover, the mass of the bottom quark is relevant in 
the light sneutrino case because if mp lT < m^, annihilation into bb is kinematically forbidden. 

3.4 Experimental constraints entering the likelihood 

We confront our model with the observables listed in Table 3. Below we comment on the 
various constraints. 

3.4.1 Relic density of sneutrinos 

We assume the standard freeze-out picture for computing the sneutrino relic abundance. The 
main annihilation channels for mixed sneutrino dark matter are i) v\i>\ — > vv (y\v\ — > vv) 
through neutralino t-channel exchange, ii) v\0\ — > ff through s-channel Z exchange, and 
Hi) v\v\ — > bb through s-channel exchange of a light Higgs. Moreover, if the v\ is heavy 
enough, it can also annihilate into W + W~ (dominant), ZZ or tt. Note that for the heavy 
LSP the annihilation into neutrino pairs is always much suppressed while the annihilation 
into other channels can be enhanced by the heavy scalar Higgs resonance. 
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i 


Nuisance parameter 
\ 


Experimental result 
A* 


Likelihood function 
Ci 


1 


m u /m d 


0.553 ± 0.043 [64] 


Gaussian 


2 


m s /m d 


18.9 ± 0.8 [64] 


Gaussian 


3 




44 ± 5 MeV [63] 


Gaussian 


4 


cr s 


21 ±7 MeV [63] 


Gaussian 


5 


P~DM 


0.3 ±0.1 GeV/cm 3 [65] 
0.43 ±0.15 GeV/cm 3 [66] 
=> 0.34 ± 0.09 GeV/cm 3 


Weighted Gaussian average 


6 


vo 


242 ± 12 km/s [67] 
239 ± 11 km/s [68] 
221 ± 18 km/s [69] 
225 ± 29 km/s [70] 
=> 236 ± 8 km/s 


Weighted Gaussian average 


7 


Vesc 


550 ± 35 km/s [71] 


Gaussian 


8 


m t 


173.3 ± 1.1 GeV [72] 


Gaussian 


9 


m b {m b ) 


4.19^* GeV [73] 


Two-sided Gaussian 


10 


a s (M z ) 


0.1184 ±0.0007 [73] 


Gaussian 



Table 2. Nuisance parameters in the scan. The values of the astrophysical parameters are taken 
from Ref. [59]. 



The annihilation into neutrino pairs proceeds mainly through the wino component of 
the t-channel neutralino and is proportional to sin it is largest for light winos. The Z 
exchange is also proportional to sin Of). The light Higgs exchange, on the other hand, is 
proportional to (A^ sin$£>)^. The dependence of 17/i 2 on the sneutrino mass and mixing angle 
has been analyzed in [40, 44]. 

We assume a 10% theory uncertainty on Qh 2 , mostly to account for unknown higher- 
order effects. In the light DM cases, one also has to worry about the change in the number 
of effective degrees of freedom in the early Universe, g e g, especially when moM ~ 20Tqcd- 
While we do take into account the change of g e g in the calculation of the relic density, the 
uncertainty related to it is not accounted for separately. Rather, we assume that it falls 
within the overall 10% theory uncertainty. (We discuss the issue of <? e ff in more detail in 
Appendix A.) 

The same annihilation channels will be relevant for indirect DM detection experiments, 
looking for gamma-rays (Fermi-LAT, H.E.S.S.), charged particles (positrons, antiprotons; 
PAMELA, Fermi-LAT, AMS) or neutrinos (Super-Kamiokande, IceCube, ANTARES), that 
could be produced by annihilation of dark matter, especially in high density regions, see 
Section 3.4.6. 

3.4.2 Direct detection limits 

The spin-independent (SI) scattering of v\ on nucleons occurs through Z or Higgs exchange. 
The Z exchange is again suppressed by the sneutrino mixing angle, while the Higgs exchange 
is enhanced by the A„ term. A peculiarity of the Z-exchange contribution is that the proton 
cross section is much smaller than the neutron one, with the ratio of amplitudes f p /f n = 
(1 — 4 sin 2 9\y). The Higgs contribution on the other hand, which becomes dominant for large 
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i 


Observable 

m 


Experimental result 
A 


Likelihood function 


1 




0.1123 ±0.0035 [46] 
(augmented by 10% theory uncertainty) 


Gaussian 


2 


UN 


(moMi &n) constraints from 
XENON10 [74], XENON100 [75], 
CDMS [76] and CoGeNT [77] 


£ 2 = e -XD D / 2 


3 


Ar z 


< 2 MeV (95% CL) [78] 


£ 3 = F( M3 ,2MeV) 


4 


Higgs mass 
limits 


from 

HiggsBounds 3.6.1beta [57, 58] 


£4 = 1 if allowed 
£ 4 = 10" 9 if not 


5 


m xt 


> 100 GeV [79] 


£5 = 1 if allowed 
£ 5 = IO- 9 if not 


6 




> 100 GeV [80] 


Cq = 1 if allowed 
£ 6 = l(r 9 if not 


7 


m fl 


> 85 GeV [80] 


£7 = 1 if allowed 
£ 7 = l(T 9 if not 


8 


nig 


> 750, 1000 GeV [81, 82] 
or none 


not included 
(a posteriori cut) 


9 


B(b -»• sj) 


(3.55 ± 0.34) x 10" 4 [83, 84] 


Gaussian 


10 




< 1.26 x 10" s (95% CL) [85, 86] 


F(/i 10 ,1.26 x 10-«) 


11 


Aa M 


(26.1 ± 12.8) x 10" 1U [87-89] 


Gaussian 



Table 3. Experimental constraints used to construct the likelihood. Where relevant, experimental 
and theoretical uncertainties are added in quadrature; in particular for Qh 2 we assume an overall 
uncertainty of (0.0035 2 + 0.01123 2 ) 1 / 2 = 0.0118. 



values of A„, is roughly the same for protons and neutrons. The total SI cross section on a 
nucleus N is obtained after averaging over the v\N and v^N cross sections, where we assume 
equal numbers of sneutrinos and anti-sneutrinos. We note that the interference between the 
Z and h° exchange diagrams has opposite sign for D\N and v*N, leading to an asymmetry 
in sneutrino and anti-sneutrino scattering if both Z and Higgs exchange are important. All 
these effects are taken into account when we compute the normalized scattering cross section 

UN- 

(Zf p + (A-Z)f n ) 2 
a N - — — 2 , (3.5) 

where fi x is the sneutrino-nucleon reduced mass, Z is the atomic number and A the mass 
number. This cross section can be directly compared to the experimental limits on cr^ 1 , which 
are extracted from the observed limits on the LSP-nucleus scattering cross section assuming 
fp = fn- 

We consider the limits coming from various direct detection experiments. In particular, 
we take into account the light dark matter results from XENON10 [74] and CDMS [76], as 
well as the latest XENON100 [75] and CoGeNT [77] results. Thus, we are using the best limits 
from both low and high mass regions, with Xenon (XENON10/100) and Germanium (CDMS 
and CoGeNT) detectors. We include the data from these experiments using a private code 
based on Refs. [90-92], where further details on the analysis can be found. For XENON100 
we adopt the best-fit light-yield efficiency L e g curve from [75]. Especially for the low DM 
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mass region, the energy resolution close to the threshold is important. We take into account 
the energy resolution due to Poisson fluctuations of the number of single electrons. The 
XENON10 analysis is based on the so-called S2 ionization signal which allows to go to a 
rather low threshold. In this case we follow the conservative approach of [74] and impose a 
sharp cut-off of the efficiency below the threshold, which excludes the possibility of upward 
fluctuations of a signal from below the threshold. Our analysis tries to approximate as close 
as possible the one performed in [74]. From CDMS we use results from an analysis of Ge 
data with a threshold as low as 2 keV [76]. We use the binned data from Fig. 1 of [76] and 
build a x 2 ) where we only take into account bins where the predicted rate is larger than the 
observed data. This ensures that only an upper bound is set on the cross section. We proceed 
for CoGeNT in a similar way. We ignore the possibility that hints for an annual modulation 
in CoGeNT are due to DM (see also [93]), and use a similar x 2 method as for CDMS to set 
an upper bound on the scattering cross section. The code allows for a consistent variation of 
the astrophysical parameters Vq, v esc and /Odm for all considered experiments. 

The information from DD is included in the Bayesian analysis in the following way. 
For XENON10 and XENON100 data, we apply the so-called maximum-gap method [94] to 
calculate an upper bound on the scattering cross section for a given mass. The probability 
returned by the maximum-gap method as a function of the model parameters as well as 
astrophysical parameters (appropriately normalized) is considered as the likelihood function 
which then is converted into the posterior PDF within the Bayesian analysis. This is an 
approximation to a pure Bayesian treatment with the advantage that it allows us to use the 
maximum-gap method, which offers a conservative way to set a limit in the presence of an 
unknown background. Since the shape of the expected background distribution is neither 
provided for XENON10 nor XENON100, it is not possible to construct a "true" likelihood 
from the data and we stick to the above mentioned approximation based on the maximum- 
gap method. 4 For CDMS and CoGeNT, the likelihood is obtained from the individual x 2 
functions as C oc exp(— x 2 /2). The method to construct the x 2 described in the previous 
paragraph amounts to introducing the unknown background in each bin i as a nuisance pa- 
rameter bi which is allowed to vary by maximizing the likelihood function under the condition 
hi > 0. Again this is an approximation to a pure Bayesian approach (in which the posterior 
PDF would be integrated over the nuisance parameters), which suffices for our purpose. 

3.4.3 Z invisible width 

A light sneutrino with ma < Mz /2 will contribute to the invisible width of the Z boson, well 
measured at LEP [78], thus putting a constraint on the sneutrino mixing: 



where T u = 166 MeV is the partial width into one neutrino flavor. For one light sneutrino with 
mp-L = 5 (20) GeV, this leads only to a mild constraint on the mixing angle of smO^ < 0.39 
(0.43). For rriu lT = 4 GeV, mp^ = = 5 GeV and assuming a common mixing angle, 

this constraint becomes stricter: sm9„ < 0.3. 

4 In [95] XENON100 data has been implemented in a Bayesian study by constructing a likelihood function 
from the Poisson distribution based on the total number of expected signal and background events. We have 
checked that such a procedure leads to similar results as our approach based on the maximum-gap method. 




(3.6) 
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On the other hand, a minimum amount of mixing is needed for light Pi's to achieve 
large enough annihilation cross section. In [40] we found sinOa > 0.12 for LSP masses above 
the 6-threshold, where annihilation into bb through Z or h can contribute significantly, and 
s'mOo > 0.25 for niy x < m^. Therefore, for light sneutrinos, the mixing angle should be 
not far from the limit imposed by the Z invisible width. Such a large mixing is however in 
conflict with DD limits unless mo l < 7 GeV. For sneutrino LSPs with masses of, roughly, 
7-40 GeV, the DD limits constrain sin 0j) to be smaller than about 0.05~ 0.07, which makes it 
impossible to achieve low enough Qh 2 . For heavier masses, one needs near the Higgs pole 
or above the W + W~ threshold to satisfy both DD and relic density constraints. This was 
also discussed in [44]. As mentioned, this splits our parameter space into two distinct regions 
where the Markov Chains converge, one with mo 1 < 7 GeV and one with > Mz/2 (more 
precisely, > 50 GeV). 

3.4.4 Higgs and SUSY mass limits 

In the MCMC sampling, we impose chargino and charged slepton mass limits [79, 80] from 
LEP as listed in Table 3. We here choose conservative values because the LEP analyses 
in principle assumed a neutralino LSP, and hence the parametrization of the LEP limits in 
terms of e.g. the chargino-neutralino mass difference as implemented in micrOMEGAs does 
not apply. To evaluate Higgs mass constraints based on LEP, Tevatron and LHC data, we 
use HiggsBounds 3 . 6 . lbeta. (The latest CMS limit on A/H — > tt [96] is also included 
via HiggsBounds.) Here note that for large sneutrino mixing, which as detailed above is 
necessary for light mo 1 , the light Higgs mass receives an important negative correction from 
the sneutrino loop, which is proportional to {A^ 4 / '(m~ 2 — m~^) 2 . Thus the lower limit on m h o 
also somewhat constrains the sneutrino sector. In order to take into account the theoretical 
uncertainty in m h o, we smear the Higgs mass computed with SuSpect by a Gaussian with 
a width of 1.5 GeV before feeding it to HiggsBounds. In the light sneutrino case, the Higgs 
decays into sneutrinos are always kinematically allowed, and they are enhanced as A„; as a 
result the h° decays almost completely invisibly in this case. (In the heavy sneutrino case, 
only a small fraction of the points have < m^o/2.) The Higgs decays into sneutrinos are 
properly taken into account in our HiggsBounds interface. 

An important point of our study is how SUSY mass limits from the 2011 LHC searches 
affect the sneutrino DM scenarios. Here note that squarks and gluinos undergo the same 
cascade decays into charginos and neutralinos as in the conventional MSSM. Since we assume 
gaugino mass unification, the gluino and squark mass limits derived in the CMSSM limits 
from jets+i?™ lss searches apply to good approximation. We have checked several v\ LSP 
benchmark points and found m g > 750 GeV for niq ~ 2 TeV based on a fast simulation 
of the ATLAS 0-lepton analysis for 1 fb _1 [81]. This is in very good agreement with the 
corresponding gluino mass limit in the CMSSM for large niQ. For 5 fb _1 of data, this limit 
should improve to m g > 1 TeV. 

A word of caution is in order however. For niq 3> nig we expect g — > qqxi 2 an d 
g — > qq')ct as hi t ne MSSM with a neutralino LSP. In our model, the Xi 2 decay further into 
the v\ LSP; if this decay is direct, x± : 2 ~~ v ^li it is completely invisible. Indeed, the invisible 
Xi 2 decays often have close to 100% branching ratio. We do not expect however that this has 
a large effect on the exclusion limits. The situation is different for chargino decays. In most 
cases, the xf decays dominantly into a sneutrino and a charged lepton (e, u or t, depending 
on the sneutrino flavor). This can lead to a much larger rate of single lepton or dilepton 
events. As a consequence, we expect the limits from 0-lepton jets+I?™ lss searches to weaken, 
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while single lepton or dilepton +E™ 1SS searches should become more effective than in the 
CMSSM. Overall, assuming gaugino mass unification, the gluino mass limit should remain 
comparable to the limit derived in the CMSSM. 

A detailed analysis of the SUSY mass limits in the sneutrino DM model is left for a 
separate work. In the present paper, we are interested in the effect of the LHC pushing the 
gluino mass limit to nig > 750 GeV or m g > 1000 GeV, see above. In order to illustrate 
this effect without having to run the MCMC several times (which would have been too CPU 
intensive), we apply the gluino mass constraint a posteriori. As we will see, it is only relevant 
for the light sneutrino case. 

3.4.5 Low-energy observables 

Further important constraints on the model come from flavor physics and from the muon 
anomalous magnetic moment. Regarding flavor physics constraints, we use the HFAG aver- 
age value of B(b -> sj) = (3.55 ± 0.24 ± 0.09) x 10~ 4 [83] with a theoretical uncertainty of 
0.23 x 10~ 4 [84]. Moreover, we use the combined LHCb and CMS limit on B(B S -> [85], 
augmented by a 20% theory uncertainty (mainly due to fs a ) as suggested in [86] . After com- 
pletion of the MCMC runs, a new limit of B(B S -> n + fT) < 4.5 x 10~ 9 (95% CL) [97] became 
available. We impose this new limit a posteriori, again assuming 20% theory uncertainty, 
but the effect of this on the posterior distributions is marginal. 5 

Regarding the supersymmetric contribution to the anomalous magnetic moment of the 
muon, Aa^ USY , we implement the 1-loop calculation taking into account the mixing between 
RH and LH v^. Then we compare this value to Aa^ = a^ xp — a^ M , where for a^ xp we take the 
experimental value reported by the E821 experiment [88], and for a^ M we take the result of 
Ref. [87] (note however the slightly lower a^ M reported in [98]). Guided by [89] and because 
of our ignorance of the 2-loop effects involving mixed sneutrinos, we assume a conservative 
theoretical uncertainty of 10 x 10~ 10 . This brings us to Aa^ USY = (26.1 ± 12.8) x 10~ 10 in 
Table 3. 

3.4.6 Indirect detection of photons and antiprotons 

Dwarf Spheroidal galaxies (dSphs) in the Milky Way provide a good probe of DM through 
the observation of gamma-rays. Although the photon signal is weaker than from the Galactic 
center, the signal-to-noise ratio is more favorable since dSphs are DM dominated and the 
background from astrophysical sources is small. From measurements of the gamma-rays from 
ten different dSphs [99], the Fermi-LAT collaboration has extracted an upper limit on the 
DM annihilation cross section in three different channels: W + W~ , bb, and t + t~. For this 
one assumes a NFW [100] dark matter profile. For DM lighter than 40 GeV, both the bb 
and t + t~ channels have the sensitivity to probe the canonical DM annihilation cross section, 
av > 3 x 10~ 26 cm 3 /s. We will not use these constraints in the fit but rather compare our 
predictions for the annihilation cross section in different channels with the limit provided 
by Fermi-LAT . We will see in the next section that this measurement constrains sneutrino 
DM in only a few scenarios for three reasons. First, for light sneutrinos we have a sizable 
v\ {y\) pair annihilation into vv (vi>), which clearly cannot lead to a photon signal. Second, 
Fermi-LAT has not published results for DM particles lighter than 5 GeV, where the bulk 
of our light DM sample that survives direct detection constraints lies. Third, Fermi-LAT's 

Effectively, we impose B(B S — > fi + fx~) < 5.4 x 10~ 9 as a hard cut, but we have checked that this makes 
no difference as compared to reweighing the likelihood according to eq. (3.4). 
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sensitivity is still one order of magnitude above the canonical cross section for DM masses at 
the electroweak scale or above. 

Annihilation of DM in the Milky Way will also, after hadronisation of the decay prod- 
ucts of SM particles, lead to antiprotons. This antiproton flux has been measured by 
PAMELA [101] and fits rather well the astrophysics background [102]. There is however 
a large uncertainty in the background at low energies (below a few GeV) due to solar mod- 
ulation effects that are not well known. Furthermore the antiprotons — as well as any other 
charged particle — propagate through the Galactic halo and their energy spectrum at the 
Earth differs from the one produced at the source. The propagation model introduces addi- 
tional model dependence in the prediction of the antiproton flux from DM annihilation. As 
for photons above, we will not use the antiproton flux as a constraint in the fit, but com- 
pare our predictions for different propagation model parameters with the measurements of 
PAMELA. We will see that the largest flux, and the largest deviation from the background, 
are observed at low energies when the sneutrino DM has a mass of a few GeV, thus leading 
to an excess of events for some values of the propagation parameters. 

Finally, a comment is in order regarding annihilation into neutrinos. Indeed, neutrino 
telescopes (Super-Kamiokande, IceCube, ANTARES) may probe sneutrino DM annihilation 
into neutrinos, e.g. from the Galactic Center or from accretion in the Sun. The neutrino 
flux from annihilation of DM captured by the Sun is determined by the cross section for 
sneutrino scattering on nucleons discussed in [40] and Section 3.4.2. We do not include a 
possible neutrino signal in this analysis but leave it for a future study. 

4 Results 

Let us now present the results of this analysis. As mentioned, for each of the three scenarios 
which we study, we run 8 Markov chains with 10 6 iterations each. The distributions of the 
points in these chains map the likelihood of the parameter space. We hence present our 
results in terms of posterior probability distributions shown in the form of histograms (1- 
dimensional distributions) with 100 bins and of contour graphs (2-dimensional distributions) 
with 100 x 100 bins. Results based on alternative (logarithmic) priors in the sneutrino sector 
can be found in Appendix C. 

4.1 Light sneutrino DM with mass below 10 GeV 

We begin with the case of light sneutrinos that was previously studied by some of us in [40] . 
Figure 1 shows the 1-dimensional (ID) marginalized posterior PDFs of various interesting 
quantities, in particular sneutrino masses and mixing angles, A terms, squarks, gluino and 
Higgs masses, etc. The blue histograms are the posterior PDFs taking into account con- 
straints 1-7 and 9-11 of Table 3, while the black (red) lines show the posterior distributions 
after requiring in addition that the gluino be heavier than 750 (1000) GeV. Note that a lower 
bound on the gluino mass not only cuts the peak of the gluino distribution but also leads to 
a lower bound on the chargino and neutralino masses, since 6m^o ~ 3m^+ « nig. (We do 
not show the m^o, m^o, m~± posterior probabilities in Fig. 1, because they follow completely 
the rrig distribution.) 

As can be seen, the DD limits, in particular from XENON10, require the sneutrino LSP 
to be lighter than about 7 GeV, with the distribution peaking around 4 GeV. (The shoulder 
at 4.5-5 GeV is due to the onset of the bb annihilation channel.) For LSP masses below 
4 GeV, the DD limits are not important. Indeed the largest cross section, obtained with the 
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Figure 1. Posterior PDFs in ID for the light sneutrino case. Specific values for best fit and quasi- 
mean points as well as the 68% and 95% BC intervals are given in Appendix B. 
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Figure 2. Posterior PDF of sin 6^ versus m^ lT for the light sneutrino case. The black and grey 
lines show the 68% and 95% BCRs before gluino mass limits from the LHC. The red and blue regions 
are the 68% and 95% BCRs requiring rrig > 1 TeV. The green star marks the bin with the highest 
posterior probability after the gluino mass limit, while the yellow diamond marks the mean of the 2D 
PDF. The grey star/diamond are the highest posterior and mean points before imposing the gluino 
mass limit. 

maximum value of sin 9„ T allowed by the Z invisible width, is below the current experimental 
limits [40]. The gluino mass bound from the LHC disfavors very light sneutrinos of about 
1-3 GeV, because the tqi'i — > vv and v\v\ — > vv annihilation channels get suppressed (recall 
that we assume GUT relations between gaugino masses). This means one needs to rely on 
annihilation through Z or Higgs exchange, as is reflected in the change of the sin Qy T and A„ T 
probability densities in Fig. 1. 

The other distributions are basically unaffected by the gluino mass cut, the exceptions 
being At and m h o. Larger values of At are preferred for nig > 1 TeV, because it is needed 
to compensate the negative loop correction to m^o from the larger Ad t in order to still 
have rriyfi > 114 GeV. Regarding m^o, the distribution is shifted towards the lower limit of 
114 GeV because of this negative loop correction. Finally, we note that the light Higgs decays 
practically 100% invisibly into sneutrinos. Therefore, should the excess in events hinting at 
a Higgs near 125 GeV be confirmed, the light sneutrino DM scenario would be ruled out. 

Regarding the supersymmetric contribution to Ao^, shown in the bottom right panel in 
Fig. 1, this is peaked towards small values. Nevertheless, the probability of falling within the 
experimental la band is sizable, p(Ao^ = (26.1 ± 12.8) x 10~ 10 ) = 31%. The larger values of 
Aa^ are obtained when there is a large contribution from the sneutrino exchange diagram. 

Our expectations regarding the relation between mass and mixing angle are confirmed 
in Fig. 2, which shows the 2-dimensional (2D) posterior PDF of sin 6^ versus ma lT . To be 
more precise, what is shown are the 68% and 95% Bayesian credible regions (BCRs) before 
and after a gluino mass cut of m g > 1 TeV. As can be seen, the region of m„ lT ~ 1-3 GeV, 
which requires sin#£ T ~ 0.3-0.4 to be consistent with WMAP, gets completely disfavored by 
a heavy gluino. 6 

In Fig. 3, we show the influence of the gluino mass limit on the predicted DD cross 
section for Xenon (we display the Xenon cross section to directly compare with the best limit 
which comes from XENON10). Imposing nig > 1 TeV has quite a striking effect, limiting 

6 To be more precise, it gets disfavored by a heavy wino, since m§ > 1 TeV implies m^o > 300 GeV in our 
model. 
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Figure 3. On the left, 2D posterior PDF of axe versus mo lT before and after imposing rag > 1 TeV; 
see the caption of Fig. 2 for the meaning of colors and symbols. On the right, correlation between 
crxe and gluino mass; the red and blue areas are the 68% and 95% BCRs. 




Figure 4. 68% and 95% BCRs of v esc versus m i>lT (left) and of pdm versus m^ lT (right). The black 
(grey) contours are the 68% (95%) BCRs without gluino mass cut, while the red (blue) areas are the 
68% (95%) BCRs for rrig > 1 TeV. The dashed lines mark the la experimental bounds for v csc and 
Pdm- 



o"Xe to a small region just below the current limit. We recall that XENON10 only constrains 
the mass range above ~ 4 GeV; for lower D\ masses, the DD cross section is constrained 
from above by the Z invisible width. We also note that there is a lower limit on the DD 
cross section [40], so that if a lower threshold can be achieved to probe masses below 4 GeV, 
in principle the light sneutrino DM case can be tested completely. (For mpj ~ 4-6 GeV, 
an improvement of the current sensitivity by about a factor 3 is sufficient to cover the 95% 
region, while an improvement by an order of magnitude will completely cover this mass 
range.) 

The influence of the nuisance parameters is also interesting. For example, a low local 
DM density can bring points with high DD cross section in agreement with the XENON10 
limits. Likewise, a small mixing angle at sneutrino masses around 4 GeV allows for higher 
Pdm, because the DD cross section is low. Analogous arguments hold for vo and v esc , since for 
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Figure 5. 68% and 95% BCRs for av versus sneutrino mass in various channels. Color code as in 
the previous figures. The dashed lines correspond to the Fermi-LAT limit [99], where for cc we have 
used the same value as for bb. Note that for mp 1 < m^, the cross section is zero, however to display 
this region we have arbitrarily set it to <jv b i = 10~ 31 cm 3 /s. 



light DM one is very sensitive to the tail of the velocity distribution. The effect is illustrated 
in Fig. 4. 

The MCMC approach also permits us to make predictions for the annihilation cross 
section of light sneutrino dark matter into different final states, relevant for indirect DM 
searches, see Fig. 5. When > m&, the dominant DM annihilation channels are into vv 
or bb pairs. The latter will lead to a large photon flux — in fact the partial cross section into 
bb is always in the region constrained by Fermi-LAT when mo 1 > 5.2 GeV. 

For lighter DM, the charged fermions final states giving photons are cc and t + t~. Here 
note that for a given LSP mass, imposing the lower limit on the gluino mass selects the upper 
range for both av C c and av T + T - while having only a mild effect on crv^. In particular the cc 
channel typically has a large cross section of av C c > 10 -26 cm 3 /s throughout the 95% BCR 
when rrig > 1 TeV. This could hence be probed if the Fermi-LAT search was extended to a 
lower mass range. 

Regarding annihilation into neutrinos, as mentioned earlier, the gluino mass limit strongly 
constrains scenario where annihilation into neutrino pairs is dominant, leading to an upper 
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limit of av uu + avpp < 1 x 10 -26 cm 3 /s, see the bottom-right panel in Fig. 5. A discussion 
of the neutrino signal for light sneutrino DM can be found in [40]. As mentioned, we leave a 
more detailed analysis of neutrinos from the Sun for a future work. 

Dark matter annihilation in our galaxy can also lead to antiprotons. To illustrate the 
impact of the antiproton measurements on the parameter space of the model, we have com- 
puted the antiproton flux for some sample points and compared those to the flux measured 
by PAMELA [101]. To compute this flux we have used the semi-analytical two-zone propa- 
gation model of [103, 104] with two sets of propagation parameters called MIN and MED, 
see [54]. For the background we have used the semi-analytical formulas of [102] with a solar 
modulation of = 560 MeV, which fit well the measured spectrum of PAMELA. 

The first sample point has a DM mass of 4.8 GeV and is dominated by annihilation 
into bb with av^ = 1.1 x 10~ 26 cm 3 /s. The resulting antiproton flux is displayed as the blue 
band in Fig. 6. A large excess is expected at energies below 1 GeV for MED propagation 
parameters, corresponding to the upper edge of the blue band. With MIN propagation 
parameters however, the flux exceeds the la range only in the lowest energy bin (Ep = 
0.28 GeV). We therefore conclude that such sneutrino DM would be compatible with the 
PAMELA measurements only for a restricted choice of propagation model parameters. Here 
note that the lowest energy bins are the ones where the background is most affected by solar 
modulation effects. 

The second sample point has lighter DM, = 2.3 GeV, and annihilation into c- 
quarks dominates the hadronic channels (av C c = 1.7 x 10~ 26 cm 3 /s) although the dominant 
annihilation channel is into neutrinos. The antiproton flux is therefore expected to be both 
lower and shifted towards lower energies as compared to the previous case. We find that the 
antiproton flux again exceeds the measured spectrum by more than la only in the first energy 
bin. Such a sneutrino DM is therefore not constrained by the antiproton measurements unless 
one chooses propagation parameters that lead to large fluxes. In this respect note that we 
can of course get even larger fluxes than those displayed in Fig. 6 using the MAX set of 
propagation parameters. 
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Figure 6. Antiproton flux as a function of the kinetic energy of the antiproton for two representative 
points as described in the text. The blue (yellow) band corresponds to n%v x — 4.8 (2.3) GeV, with the 
upper curve corresponding to MED and the lower curve corresponding to MIN propagation parame- 
ters. We also display the background only (black line) and the PAMELA data for energies below 10 
GeV (red crosses). 
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4.2 Heavy sneutrino DM 



Let us now turn to the case of heavy sneutrinos. We will first discuss the heavy non- 
democratic (HND) case, where the LSP is the vx T , and then the heavy democratic (HD) 
case, where all three neutrinos are close in mass and any of them can be the LSP or co-LSP. 

The posterior PDFs in ID for the HND case are shown in Fig. 7. Here, we do not 
superimpose the distributions with nig > 750 or 1000 GeV, because the gluino automatically 
turns out heavy, with 99% probability above 1 TeV. The v\ T masses now range from 90 to 
255 (80 to 375) GeV at 68% (95%) Bayesian credibility. There is also a small region near 
rriv lT ~ 60 GeV, where the sneutrino annihilates through the light Higgs resonance; this 
region has 3% probability. 7 See Table 6 in Appendix B for more details. The V2 T mass is 
typically very heavy, above 1 TeV, and the mixing angle is required to be very small to evade 
the DD limits, cf. the discussion in Section 3.4.2. Interestingly, the mixing can be almost 
vanishing; this happens either when rriu lT ~ m h o J 2 so that the annihilation is on resonance, 
or when co-annihilation channels are important. In the first case, the Ay term must be very 
small, otherwise the annihilation cross section would be too large and Qh 2 too small. Note 
that the upper limit on the sneutrino LSP mass is determined by the range for the gluino 
mass used in the scan which in turn sets an upper bound of 500 GeV on the lightest neutralino 
and hence on the sneutrino LSP. 

The light Higgs mass is not much affected by radiative corrections from a heavy sneu- 
trino, so the posterior PDF of m^o is like in the conventional MSSM. (See the bottom row 
of Fig. 7 for Higgs-related quantities.) A light Higgs in the 123-127 GeV mass range has 
21% probability in this case. As in the MSSM, this mass range requires large mixing, see the 
distribution for Xt/Ms- S The signal strength in the gg — > h — > 77 channel relative to SM 
expectations [R ggil ) is also just like in the MSSM [105], with the highest probability being 
around Rgg-yy ~ 0.9. In this scenario, it is much more difficult to reach larger values of Aa^ 
as the sneutrino contribution is never large. We find Aa M < 8.6 x 10~ 10 at 95% BC. 

In Fig. 8, we show the 2-dimensional posterior PDF of sin#j> T versus m^ T . As can be 
seen, the mixing angle is always in the sin Qy T ~ 0.01 — 0.05 region except when ~ m^o /2 
or for a few scattered points with heavier LSP masses. The latter correspond to cases where 
the co-annihilation of pairs of NLSPs nearly degenerate with the sneutrino LSP helps to 
increase the effective annihilation cross section, so that the relic density of the sneutrino is 
in agreement with WMAP. The NLSP can be either a neutralino or a slepton. For the bulk 
of the points, however, the minimal value of the mixing increases with the sneutrino mass. 

The predictions for the SI cross section are within one order of magnitude of the XENON 
and CDMS bounds except when rrij> lT ~ and for the scattered point where coannihi- 

lation dominates, see the right panel in Fig. 8. Indeed, when the annihilation in the early 
Universe is enhanced by a resonance effect, the coupling of the LSP to the Higgs has to be 
small, hence one needs a small mixing angle. This also means that the sneutrino coupling to 
the Z is small, leading to a small SI cross section. 



7 As mentioned in Section 3.4.1, the sneutrino can also annihilate through the heavy scalar (not the pseu- 
doscalar!) Higgs resonance. We have checked that this process does occur in our chains. However, it turns 
out that it is statistically insignificant and does not single out any special region of parameter space. 

8 X t = At — (jl/ tan/5 and Mj = m^m^. In fact the distribution of A t is the only one that is significantly 
changed by requiring m h a G [123, 127] GeV, see also [105]. 
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Figure 7. Posterior PDFs in ID for the heavy non-democratic (HND) sneutrino case. 
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500 



v\ T mass (GeV) v\ T mass (GeV) 

Figure 8. Posterior PDFs in 2D of sm6p T (left) and ctxc (right) versus mi> lT for the HND case. 
The red and blue areas are the 68% and 95% BCRs, respectively. The green stars mark the highest 
posterior, while the yellow diamonds mark the mean of the PDF. 

The precise relation between the LSP mass and the Higgs mass has important conse- 
quences when we consider annihilation channels in the galaxy. In some cases, such annihi- 
lations can be strongly enhanced with respect to their values in the early Universe. This 
Breit-Wigner enhancement can occur when the annihilation proceeds through a s-channel 
exchange of a Higgs particle near resonance, the cross section is then sensitive to the thermal 
kinetic energy: at small velocities, one gets the full resonance enhancement while at v ~ c, 
one only catches the tail of the resonance [106-109]. This occurs when 1 — m^ /4m? iT <C 1, 
thus when the annihilation is primarily into bb. In the left panel in Fig. 9, a small region 
at 95% BC has a photon flux above the limit imposed by Fermi-LAT . Away from this spe- 
cial kinematical configuration, the annihilation cross section into bb is usually two orders of 
magnitude below the present limit. The dominant annihilation channel is rather into W- 
boson pairs. Even for this channel, the predictions are at least one order of magnitude below 
the Fermi-LAT limit except when m„ lT ~ 100 GeV, where the predictions are only a factor 
2-3 below the limit. The annihilation into neutrino pairs is always subdominant for heavy 
sneutrinos, with av vv + avpp < 10 -30 cm 3 /s. 

Note that even after removing the points that are excluded by Fermi-LAT in the bb 
channel, the predictions for crxe extend to small values. Indeed for these points there is no 
large enhancement of the annihilation rate in the early Universe, hence no need to have small 
couplings of the LSP to the Higgs. Therefore the predictions for the SI cross section covers 
a wide range and is not correlated with crv b i, see the bottom right plot in Fig. 9. 

We have also computed the predictions for the antiproton flux for the heavy sneutrino 
case. The largest fluxes are expected for DM masses around 100 GeV where the annihilation 
cross section can reach 3x 10~ 26 cm 3 /s. We found that with the MED propagation parameters 
the flux is barely above the background and always within the la experimental error bars. 
Note that a large flux is also expected for the few points that have a large annihilation into 
bb, these points are however already excluded by Fermi-LAT as discussed above. 

The results discussed above for the HND case also hold for the HD sneutrino case. In 
fact, most of the distributions in the HD case are practically the same as in the HND case. 
The only differences are observed for the LSP mass, and for the associated see Fig. 10. 
We note a slightly higher probability of 6% to be on the h° pole. Correspondingly, also 
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Figure 9. 2D posterior PDFs for the HND case relevant for indirect DM detection; color codes etc. as 
in Fig 8. 

small Ay and small mixing angles have somewhat higher probability than in the HND case. 
Regarding the flavor of the LSP, we find that a r-sneutrino LSP has 55% probability and is 
thus, as expected, somewhat preferred over e/n sneutrino co-LSPs (45% probability), see the 
right-most panel in Fig. 10. The fact that the v\ e —v\ T mass difference peaks within ±10 GeV 
is however just a consequence of our prior assumption for the HD case. 
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Figure 10. Posterior PDFs in ID for the heavy democratic (HD) sneutrino case. All other distribu- 
tions are practically the same as in the HND case. 
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5 Conclusion 



We performed a global MCMC analysis of a sneutrino DM model with Dirac neutrino masses 
originating from supersymmetry breaking. The main feature of this model is a mainly RH 
mixed sneutrino as the LSP, which has a large coupling to the Higgs fields through a weak- 
scale trilinear A term which is not suppressed by small Dirac-neutrino Yukawa couplings. We 
demonstrated that such a RH sneutrino can be an excellent cold dark matter candidate over 
a wide range of masses. In particular, it can be consistent with all existing constraints for 
masses around 3-6 GeV, as well as for masses of about 50-500 GeV (the upper limit coming 
from the fact that we consider gluino masses only up to 3 TeV). 

Direct detection limits in particular from XENON10 heavily constrain the low mass 
range. The DD cross section however sensitively depends on, e.g., the escape velocity in 
the light DM case. We therefore took special care to account for uncertainties arising from 
astrophysical parameters, like vq, v esc and the local DM density pdm- Moreover, we accounted 
for uncertainties from the quark contents of the nucleon, relevant for the Higgs exchange 
contribution to the DD cross section. 

Our main results are posterior probability distributions of parameters, masses, and 
derived observables — in particular the LSP mass and the direct and indirect detection cross 
sections. Assuming gaugino-mass unification, the recent LHC limits on the gluino mass 
exclude the very light sneutrino DM region below about 3 GeV, where the DD limits are not 
efficient. To be precise, requiring m g > 1 TeV leads to 2.9 < < 5.6 GeV at 95% BC. For 
heavy sneutrinos of the order of 100 GeV, the gluino is always heavy so that the LHC limits 
have no effect on sneutrino DM. 

Regarding the prospects for probing light sneutrino DM, we found that covering the 
95% BC region requires about a factor three increase in sensitivity in direct detection for 
DM masses around 5 GeV, as well as a lower threshold to be able to probe masses below 
4 GeV. Similarly, the prospects for indirect detection through photons and antiprotons are 
promising if the sensitivity of experiments can be extended to lower masses. The crucial 
test however comes form the LHC: in the light sneutrino DM scenario the Higgs decays 
dominantly invisibly into sneutrinos. Therefore if the Higgs-like excess around 125 GeV is 
confirmed, the light sneutrino model is ruled out. 

The heavy sneutrino scenario can also be probed by DD experiments, this requires an 
increase in sensitivity of roughly one order of magnitude over the current limits. Only a small 
region where the sneutrino has about half the mass of the light Higgs would remain out of 
reach in this case. Such a scenario should lead to a Higgs signal that is compatible with the 
SM. Prospects of indirect detection are more challenging for the heavy sneutrino. 

Both, light and heavy, sneutrino scenarios offer distinctive LHC SUSY phenomenology. 
In particular neutralinos (typically x? an d X2) appearing in squark and gluino cascades 
can decay invisibly into the LSP. Indeed the probability for a 90% invisible decay of the 
lightest (next-to-lightest) neutralino is about 80% (50%) for the light sneutrino scenario, 
and close to 100% (30-40%) in the heavy sneutrino scenarios, see Table 4. This implies 
that there can be up to three different invisible sparticles in an event. The dominant decay 
of charginos (with a branching fraction larger than 0.5) is into a charged lepton and the 
LSP with roughly 50% probability. The charged lepton is typically a r for the light and 
heavy non-democratic scenarios or a ej \i for the heavy democratic scenario. The cascade 
decays of squarks, (Jr — > qxi — > qv\v, q_L — > QX2 ~^ QV\v, q~L — > q'xf — > q'li'i therefore 
give different amount of missing energy as compared to the MSSM. Furthermore the cascade 
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light 


HND 


HD 


B(xi -»■ inv) > 0.9 


79% 


96% 


98% 


B(x2 -»■ mv) > 0.9 


53% 


29% 


42% 


B(xf -> ^ i>u) > 0.5 


7% 


9% 


48% 


B{xi r± P ir ) > 0.5 


46% 


46% 


10% 



Table 4. Probabilities for neutralino and chargino decays in the light, HND and HD sneutrino DM 
cases, requiring wig > 1 TeV. 



decays of gluinos, g — > x}jj w ^ a ^ so gi ye a large contribution to the jets plus missing Ex 
channel while the decay of gluino pairs via a chargino will give about the same amount of 
same-sign and opposite-sign lepton pairs. Note that the alternative dominant decay mode 
of the chargino is xf ~^ W^Xi'i m this case the mass of the invisible particle could be much 
larger than the DM mass. Probabilities for x!ii x\ an d X\ decays in the light, HND and HD 
scenarios are summarized in Table 4. 

Distinctive features of the light sneutrino scenario at the LHC were investigated in 
[45] and it was shown that distributions such as lepton and jet number as well as same- 
sign/opposite-sign dilepton rates could give a distinctive signature of a light sneutrino at the 
LHC. A detailed analysis of the LHC sensitivity in the heavy sneutrino DM model, based on 
publicly available ATLAS and/or CMS results for different signal topologies, is underway. 

Note added: After the completion of this work, ATLAS and CMS announced a ~ 
5<T discovery consistent with a SM-like Higgs boson around 125-126 GeV [111, 112]. This 
excludes the light sneutrino scenario discussed in our work, while the heavy sneutrino case 
remains viable. Moreover, new XENON100 results [113] became available, constraining the 
heavy sneutrino case for sin Qy above 0.03. For flat priors, the probability to obey the new 90% 
CL exclusion limit in the HND case is 62% (54% with logarithmic priors). Our conclusions 
remain unchanged. 
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A Discussion of T QCD / 'g eS / 'h eS 



In the standard freeze-out picture, Qh 2 is inversely proportional to the number of effective 
degrees of freedom, g e g. At the temperature where the QCD confinement occurs, around 
Tqcd ~ 300 MeV, g e ff starts to drop and Qh 2 increases. This is relevant for DM masses 
below ca. 7 GeV, where the freeze-out temperature Tf « mDfvi/20 is of the order of Tqcd- 

The uncertainty in the equation of state (actually in the effective degrees of freedom 
g c s and h e g contributing to the energy and entropy densities of the SM) at temperatures 
around Tqcd induces a non-negligible uncertainty in the calculated Clh 2 for light DM. In 
[110], Hindmarsh and Philipsen estimated this uncertainty to be at around the 15% level. 
They also provided five tables (A, B, B2, B3 and C) describing the evolution of the effective 
degrees of freedom in the early Universe using different parameters for the equation of state 
and different values for the temperature at which there is a sharp switch between quarks and 
gluons, and hadrons and their resonances. 

Taking 5000 sample points from our light sneutrino sample we compute VLh? using the 
five tables of [110] and compare it to the "default" value Odef^ 2 obtained with the default 
micrOMEGAs table. The result is shown in Fig. 11. We note that the variations in the 
computed VLh 2 (relative to the default value of micrOMEGAs) can be as large as 20%. The 
main upward fluctuation is due to table A, which corresponds to a case where hadrons are 
ignored in the confined phase. 

We do not take this into account as additional uncertainty in our analysis, as a somewhat 
larger uncertainty in VLh? for DM masses below about 5 GeV would not sensitively impact 
our results. However, we note that the situation is unsatisfactory and would merit further 
study. 




v\ T mass (GeV) 



Figure 11. Maximum variations of ilh 2 for light sneutrinos, due to different evolution of the effective 
degrees of freedom from the tables of [110] as compared to the default table in micrOMEGAs. Unlike 
the other tables, Table A ignores hadrons in the confined phase. 
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B Tables: Bayesian credible intervals for parameters and observables 

We provide in Tables 5, 6 and 7 the 68% and 95% BCIs (Bayesian credible intervals), which 
we define as highest posterior density intervals, for several parameters and observables, in the 
case of light, HND (heavy non-democratic) and HD (heavy democratic) sneutrinos. For each 
case, we also provide information on the "best fit point" (the point with highest likelihood in 
the MCMC sampling), and a so-called "quasi-mean point" which is close to the mean of our 
parameters. Due to correlations between parameters in order to respect the constraints, as 
well as asymmetric and multimodal distributions, the mean point itself is very unlikely. Thus 
we pick in our samples the closest point to the mean with a good likelihood as an example 
of a typical point. The SLHA files for these points are available as supplementary material 
on the arXiv. 
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Table 5. 68% and 95% BCIs and values of two example points for various parameters and observables 
in the light sneutrino case with rrig > 1 TeV. 
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Table 6. 68% and 95% BCIs and values of two example points for various parameters and observables 
in the heavy non-democratic sncutrino case. 
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Table 7. 68% and 95% BCIs and values of two example points for various parameters and observables 
in the heavy democratic sneutrino case. 
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C Logarithmic priors for the sneutrino parameters 

An alternative to the uniform prior is to use a logarithmic prior, in which all orders of 
magnitude are equally likely. In fact, one may argue that the log prior is the least informative, 
i.e. the more objective, prior associated with a dimensionful quantity. This can be obtained 
from the Jeffreys prior based on Fisher information [114]. 

To probe the impact of the log priors on our sampling, we again run 8 chains with 
10 6 iterations each, in the light sneutrino and in the heavy non-democratic sneutrino case, 
but assuming logarithmic priors on the sneutrino mass parameters mj>.. For all the other 
parameters, including the mixing angle sin 9^, we assume an uniform prior as before. We 
stress that the sneutrino ^4-terms, A^, are derived from the sneutrino masses (see eq. (2.6)) 
and are thus sensitive to the change of prior. 

In Fig. 12, we show the 2-dimensional posterior PDFs of sin^ versus mv lT and of 
o"xe versus mp lT in the light case. As expected, light sneutrino masses are favored, and as 
a consequence larger mixing angles are favored (because of m^ 2T being lighter on average, 
which implies a smaller Ap T ). With the post-LHC gluino mass limit, however, the PDFs for 
log prior and uniform prior are quite similar, cf. Figs. 2 and 3. 

Figure 13 is the same as Fig. 12 but for the HND case. The main change as compared 
to the uniform prior case, see Fig. 8, is again the preference for lighter sneutrino masses, 
roughly m$ lT < 250 GeV instead of mp lT < 375 GeV at 95% BC. Furthermore, the light 
Higgs resonance region has a probability of 33% in the log prior case, as compared to 3% in 
the case of uniform prior. 

We note the extension of the 95% BCR to low mixing angles and scattering cross sections 
around 100 GeV: it results from co-annihilation, mainly with the NLSP sneutrinos (also being 
lighter due to the log prior), but also with a light neutralino or stau. The constraint on the 
gluino mass from the LHC remains largely irrelevant in the HND case: p(rrig > 1 TeV) = 94% 
instead of 99% with uniform priors. 

In summary, imposing log priors in the sneutrino sector does not lead to dramatic 
changes, however it highlights the light Higgs resonance region and the various co-annihilation 
possibilities in the heavy non-democratic case. These specific cases could lead to a sizable 




v\ T mass (GeV) u\ T mass (GeV) 



Figure 12. Posterior PDFs in 2D of sin0p T (left) and axe (right) versus m^ lT for the LD case using 
logarithmic priors. The red and blue areas are the 68% and 95% BCRs, respectively. The green stars 
mark the highest posterior, while the yellow diamonds mark the mean of the PDF. 
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Figure 13. Posterior PDFs in 2D of sm.9o T (left) and oxc (right) versus rnp lT for the HND case using 
logarithmic priors. The red and blue areas are the 68% and 95% BCRs, respectively. The green stars 
mark the highest posterior, while the yellow diamonds mark the mean of the PDF. 

decrease of the scattering and annihilation cross sections thus making it difficult to test the 
model with future direct and indirect detection experiments. 
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